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Abstract 

The dynamics of Vlasov kinetic moments is shown to be Lie-Poisson on the dual Lie algebra 
of symmetric contravariant tensor fields. The corresponding Lie bracket is identified with 
the symmetric Schouten bracket and the moment Lie algebra is related with a bundle of 
bosonic Fock spaces, where creation and annihilation operators are used to construct the 
cold plasma closure. Kinetic moments arc also shown to define a momentum map, which is 
infinitesimally equivariant. This momentum map is the dual of a Lie algebra homomorphism, 
defined through the Schouten bracket. Finally the moment Lie-Poisson bracket is extended 
to anisotropic interactions. 
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1 Geometric methods in kinetic theory 



Liouville equation. The equations of kinetic theory in non-equilibrium statistical mechanics 
always describe the evolution of a probability distribution in phase space. The most fundamental 
kinetic equation is the Liouville equation describing the evolution of an iV-particle distribution 
p(z\, . . . ,zn), where Zi = (qi,Pi) is the phase space coordinate of the i-th. particle. Given the 
conservation of particles and the Hamiltonian nature of their dynamics, the preservation of the 
symplectic (phase space) volume assures that the equation for p is written as the advection 
equation dtp + {p, H} = 0, where the Hamiltonian H depends on all particles and {•, •} is the 
canonical Poisson bracket. 

In geometric terms the Liouville equation is a Lie-Poisson equation [44j on the group of 
canonical transformations (symplectomorphisms) on the iV-particle phase space T*Q x • • • x T*Q 
(N times) [33] • The Lie algebra of the symplectic group consists of Hamiltonian vector fields 
X# G 3t can and is identified with the Poisson algebra of generating functions H E via the 
usual Lie algebra anti-homomorphism X#[/] = {f,H}. Thus the dual Lie algebra is identified 
with the space of densities (distributions) p G T* ~ Den, which undergo Lie-Poisson dynamics. 

The BBGKY hierarchy. In the search for approximate descriptions of the Liouville dynam- 
ics, it is common to introduce global quantities that integrate out the information on some of 
the particles (usually most of them). In particular one defines a n-particle distribution as 



These quantities are called BBGKY moments and their equations constitute a hierarchy of 
equations known as BBGKY hierarchy, or Bogoliubov-Born-Green-Kirkwood-Yvon equations. 
This hierarchy is rather complicated, although Marsden, Morrison and Weinstein [33] have 
shown that this is also a Lie-Poisson system, which is defined on the Lie algebra of canonical 
transformations that are symmetric with respect to their arguments. Thus the operation of 
taking the BBGKY moments is a Poisson map. More particularly, the reason why this happens 
is that BBGKY moments are a particular kind of momentum map |44J. 

The Vlasov equation. Suitable approximations on the equation for the single particle dis- 
tribution / := fx (the first-order BBGKY moment) yield the Boltzmann equation [TJ, which we 
write simply as 



where H is now the single-particle Hamiltonian. The right hand side collects the information 
on pairwise collisions among particles. An important approximation of the Boltzmann equation 
is the Fokker-Planck equation [18| . arising from the stochastic hypothesis of brownian motion. 
However, in many contexts it is possible to neglect the effects of collisions, so that 



This is called Vlasov equation |55j and its underlying geometric structure has been widely 

investigated [38l ESI 133 [33 IB]- ^ n particular, Marsden, Weinstein and collaborators [33 EE] 
have shown that this equation possesses a Lie-Poisson structure on the group of canonical trans- 
formations of the single-particle phase space T*Q. The explicit expression of the Lie-Poisson 
bracket is 





(1) 
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where the index V will be used below to avoid confusion in the notation. 

In the general case, the Vlasov density / depends on seven variables: six phase space coordi- 
nates plus time. This indicates that the Vlasov equation is still rather complicated and it may be 
convenient to find suitable approximations. To this purpose, one introduces the moments of 
the Vlasov distribution. The use of moments in kinetic theory was introduced by Chapmann 
and Enskog [12j . who formulated their closure of the Boltzmann equation yielding the equa- 
tions of fluid mechanics. This result showed the power of using moments and the mathematical 
properties of moments have been widely investigated since that time. 



Vlasov moments I: statistical moments The investigation of the geometric properties 
of Vlasov moments mainly arose in two very different contexts, particle beam dynamics and 
shallow water equations. However it is important to distinguish between two different classes of 
moments: statistical moments and kinetic moments. Statistical moments are defined as 



M n ,n(t) := jjp n q n f(q,p,t)dqdp. 



These quantities first arose in the study of particle beam dynamics [9jll0tl40] from the observation 

1 /2 

that the beam emittance e := (.Mo,2 -M 2,0 — -M-ii) is a laboratory parameter, which is also 



invariant. In particular, Channell, Holm, Lysenko and Scovel |10l [Tl~| [28"1 |4"T] were the first to 
consider the Lie-Poisson structure of the moments, whose explicit expression is 
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This geometric framework allowed the systematic construction of symplectic moment invariants 
in [28], a question that was also pursued by Dragt and collaborators in [13]. Special truncations 
have been studied also by Scovel and Weinstein in [48] . 



Vlasov moments II: kinetic moments. Another class of moments was used by Chapman 
and Enskog [12], i.e. the kinetic moments 

A n (q,t) := J p n f(q,p,t)dp. (3) 

The following discussion will refer to these quantities as simply "moments", unless otherwise 
specified. The Hamiltonian structure of these variables was found by Kupershmidt and Manin 
[35j in the context of Benney long waves [3]. Later, Gibbons recognized that this is a Lie- 
Poisson structure inherited from the Vlasov bracket [20]. The relation between moments and 
the symplectic Lie algebra was also known to Lebedev [36], although he did not recognize the 
connection with Vlasov dynamics. The Lie-Poisson structure for the moments is also called 
Kupershmidt- Manin structure and is explicitly written as [35] 

{F, G} = — J A m+n ^i 

The main theorem regarding moments is thus the following 

Theorem 1 (Gibbons [20]) The process of taking moments of the Vlasov distribution is a 
Poisson map, which takes the Vlasov Lie-Poisson structure to the Kupershmidt- Manin bracket. 
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This argument was also used later by Gibbons, Holm and Kupershmidt [21} 122] in obtaining the 
moment equations of chromohydro dynamics, the continuum theory of a Yang-Mills quark-gluon 
plasma. 

The Lie-Poisson structure of moment dynamics remains an interesting research topic. In- 
deed, there are many open questions, some of which are answered in the present paper. For 
example, Poisson maps relating different Lie-Poisson structures are often understood in terms of 
momentum maps: What is the corresponding momentum map for kinetic moments? And also, 
the moment Lie-Poisson bracket defines a Lie algebra: What is its geometric interpretation? 
The last question is especially important in higher dimensions. Kupershmidt proposed [33|. [34] 
a multi- index treatment, so that moments may be defined as 

A a (d,t) : = y P CT /(q,p,t)d 3 p, (5) 

where a is a multi-index such that p a := p^ 1 pZ 2 P3 3 ■ The corresponding Lie-Poisson structure 
is found to be 

<*,P 3 

where lj is a multi-index such that (lj)k = $jk ( an zero entries, except the j'-th, which is equal 
to unity). However, this notation does not address the geometric interpretation of the moments. 
Are they differential forms, special tensors, or something else? 

Plan. Section [2] of this paper shows that moments are symmetric tensor densities whose Lie 
bracket is defined by the Schouten concomitant [50] . Section [3] shows that moments define a 
momentum map and investigates the properties of this map. The moment Lie algebra is then 
shown in Section U] to possess a Fock-space structure as for Bose-Einstein statistics, and its 
creation and annihilation operators are used to construct the cold-plasma closure. In Section 
[5] the moment treatment is extended to include anisotropic interactions in accord with [21[ [22] . 
Conclusions and open questions are summarized briefly in Sectional 



5F d 5G 



<7i 



6A a dqi 8 A, 



Pi 



5G d 5F 



5A P dqi 5A a 



d 3 q, 



(6) 



2 Lie-Poisson structure of kinetic moments 

The Vlasov kinetic moments have been shown [20J to possess the Kupershmidt-Manin structure 
|35] given in This property holds in one dimension and the analogous result Q holds in 
higher dimensions, when the moments are defined as in (JSJ). The procedure to obtain the Lie 
Poisson structure for the moments begins with the Vlasov bracket ([2]) and makes use of the 
chain-rule identity 

5£ = ^S_A IL 5^ = ^ SF_ 

Sf ^ Sf SA n ^ p 6A n ' [ } 

for any functional F[f]. Direct substitution of this expression into equation ([2]) yields the Lie- 
Poisson structure for the moments (j3|). In higher dimensions, the same argument holds upon 
replacing A n with A a . 



2.1 Tensorial interpretation of the moments 

In one dimension, the moment Lie-Poisson structure identifies a Lie algebra defined by the 
square bracket in @, while the multi- index notation yields the square bracket in ([6]) in higher 
dimensions. The questions addressed in this section are: What Lie algebra is this? What is its 
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geometric nature? For these considerations, the moments are regarded as fiber integrals [U 57] 
on T*Q, the Vlasov distribution is taken to be a phase space density f dq A dp and the natural 
geometry of the cotangent bundle T*Q is used to identify the momentum coordinate p with the 
canonical one form 6 = pdq. In this way the moments are defined as 

A n = [ ® n (pdq)f(q,p,t)dqAdp = A n (q,t)® n+1 dq. (8) 
JT*Q 

Since moments belong to a dual Lie algebra, the dual variables of the moments are tensor fields of 
the form (5 n = (3 n (<l, t)® n d q - Therefore, moments in one dimension are covariant tensor-densities 
dual to the space of contravariant tensor fields, endowed with the Lie bracket in equation @. 
This conclusion was obtained earlier in [21]. However, in order to enrich this analysis we need 
to extend our arguments to higher dimensions. Following the same line of reasoning allows one 
to write moments in higher dimensions as 

A a = [ {pidq 1 ) ai {p 2 dq 2 ) a2 {p3dq 3 ) a3 f(^P,t) d 3 q A d 3 p = ^(q, t) (dq) ff d 3 q , (9) 

Jt*q 

where d 3 q is the volume element on M 3 . 

At this point we prefer to forgo the multi-index notation and to adopt the definition of the 
moments that is used in the physics literature [25]. Hence, moments become defined as tensors 
of the form 

A„(q,t) = J p n /(q,p,t) d 3 p, (10) 

where p n denotes the n-th tensor power. Inserting the basis and extending to any number N of 
spatial dimensions yields 

oo „ 

A n (q, t) = / {Pi dqT /(q, P, *) d^q A d^p 
i=l Jt ZQ 

Pil ... Pin dq^ . . . dq*» /(q, p, t) d N q A d^p 

= (A n ( q ,t)) n . n dq^®---®dq l "®d N c l . (11) 

ii,—,i n 

This definition identifies the moments with symmetric covariant tensor-densities. Also, since 
contractions of the form p n _l /3 n (q) as in ([7|) must be symmetric under permutations of indices, 
the Lie algebra variable may be identified with symmetric contravariant tensor fields. 

The duality between covariant and contravariant tensors can be seen by the following argu- 
ment. Take initially /3 n = S(3 n + A(3 n as a generic contravariant n-tensor, where S and A are 
the symmetrizing and anti-symmetrizing operators respectively: the contraction with the tensor 
power p n = Sp n yields 

p n J &(q) = 5p n -l &(q) = p n _l S/? n (q) , 

since S = S* [51j . Thus, the antisymmetric part of (3 n does not contribute to the contraction and 
Pn can be identified with its symmetric part S(3 n . The same argument may be used to identify 
the dual of symmetric covariant tensor densities with the space of symmetric contravariant 
tensors. Indeed, one has {A n ,f3 n ) = (SA n ,[3 n ) = (A n ,S(3 n ), so that again (3 n is identified with 
its symmetric part. 
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2.2 The Schouten concomitant for Vlasov kinetic moments 



At this point, since the tensor notation is used in the same fashion in any number of dimensions, 
we continue to work in terms of the symmetric tensors A n . However the Lie-Poisson structure 
for the moments is well-known only in multi-index notation |33|, [34] . What is the Lie-Poisson 
structure for the tensor moments? More particularly, the Lie algebra is now the space of sym- 
metric contravariant tensor fields. What is its Lie bracket? To answer this question, it suffices 
to insert the expression 
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into the Vlasov bracket ([2]) to obtain the following Lie-Poisson structure for the moments A n 
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, where contraction of lower and upper indexes is 



with [fi n _l V)ce, 

always assumed (sum over r in this case). The key observation now is that, since moments 
are symmetric, the antisymmetric part of the tensor expression in square brackets yields zero 
contribution. Consequently, one can rewrite the Lie Poisson bracket for the moments as 
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where S is the symmetrizer operator and the square brackets [•,•] define [50] the Schouten 
concomitant (or symmetric Schouten bracket), which is a well-known object in differential 
geometry [32 . Upon denoting the Schouten bracket as [f3 n ,a m ] = &dp n a m , one may introduce 
the infinitesimal coadjoint action ad* as its dual, 



^ (ad^ n A k , a fc _ n+ i) := ^ ( A k , ad^ a k _ n+1 ) 



(16) 



k,n=0 
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Thus, ad*j n A k yields a covariant (k — n + l)-tensor density. The moment equations may be 
obtained from (1131) and written in terms of the Hamiltonian 7i as 
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An explicit expression for ad^yl m _|_ n _i may be obtained from its definition by a direct calculation 

as 

- &dp n A m +n-i = n(fi n A V)A m+n _i +n(VJ /3 n )A m+n _i + m A m+n -i -I V/3„ , (18) 
where we use the notation 
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so that the summation is intended over all the repeated indexes. For the case of Vlasov dynamics, 
one has k = m + n — 1 and since m > 0, the contravariant indexes are all contracted, so that 
the resulting quantity is always a (symmetric) covariant tensor density. 

Thus the tensor interpretation of the moments provides a direct identification between the 
moment Lie bracket and the symmetric Schouten bracket (or concomitant). This bracket was 
known to Schouten as an invariant differential operator (substitution of covariant derivatives 
is always possible) and its relation with the polynomial algebra of the phase-space functions is 
very well known. 

2.3 History of the Schouten concomitant 

The symmetric Schouten bracket has been object of some studies in differential geometry [51 \15\ 
[53j E7j and has found applications in various aspects of quantum mechanics [HE]. Also, since 
Killing tensors on a Riemannian (or pseudo-Riemannian) manifold are symmetric contravariant 
tensors, they have been associated with the Schouten concomitant in [2j [52] and it has been 
pointed out \19\ [2] that they can be defined as the moment sub-algebra consisting of symmetric 
tensor fields commuting with the contravariant metric g -1 , under the Schouten concomitant. 

To our knowledge, the relation of the Schouten concomitant with Lie-Poisson dynamics for 
the Vlasov moments is new. However, it is perhaps not entirely unexpected. In fact, the Lie- 
Poisson bracket functional (I13p was already considered by Kirillov [31] in the context of invariant 
differential operators. This is not surprising, since Kirillov was the first to introduce the Lie- 
Poisson structure for Hamiltonian systems, also known as the Kostant-Kirillov structure. In 
particular, he noticed how this bracket functional can generate what we have called the coadjoint 
operator (adp Ak), "which, apparently, has so far not been considered", he remarked in 1982. 
(The Kupershmidt-Manin operator was known since 1977, although its association with the 
Schouten concomitant had apparently not been noticed.) What Kirillov considered were the 
cases h = 1, which is the Lie derivative, and h = k, which is often called Lagrangian Schouten 
concomitant. 

2.4 Specializations of the Schouten concomitant 

As one may verify directly, the Schuouten concomitant ad^ai among two vector fields (3\ and 
ctx is equal to the Jacobi-Lie bracket [,9i,ai]. a = ad^ai, which is equal, in turn, [H] to the 
Lie derivative £p 1 cti = \j3\, ai],^. One also obtains ad^a^ = £p 1 a n and analogously the dual 
operation ad* is (minus) the Lie derivation ad^^4 n = — £p 1 A n . Moreover, it is interesting to 
notice the chain of relations, 

(ad^ n A n , ai) = - (A n , ad Ql (5 n ) = - {A n , £ ai (5 n ) =: (A n o n , oti) , 

where (•, •) is the L 2 pairing and the diamond operation o denotes (minus) the dual action of 
the Lie derivative, following the work of Holm, Marsden and Ratiu [27]. Thus, one obtains the 
following expression for the Lagrangian Schouten concomitant: 

a d/3n A n = n o A n . 

Consequently, the equation for the first-order moment A\ (the fluid momentum) may be written 
from (TTT]) as 

-^ = 2_^Pn<>A n . 

n 

In fluid dynamics, the first two terms in the sum represent the fluid transport (n = 1) and the 
pressure pressure contribution (n = 0) for fluid motion [27] . 
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2.5 The graded structure of the moment algebra 



The Lie algebra structure generated by the Schouten concomitant is that of a graded Lie algebra 
= ©jfli- The filtration [g n ,g m ] = flm+n-i allows the subspace gi = X of vector fields as 
a particular sub-algebra, whose underlying Lie group consists of the diffeomorphisms on the 
configuration manifold [23]. The largest possible sub-algebra corresponds physically [23] to the 
barotropic fluid motion and it is given by gi ® go, where go consists of scalar functions. 

2.6 Momentum map considerations 

The Schouten concomitant establishes an isomorphism g — > T between the moment algebra g 
and the polynomials in the p's, which are embedded in the space T of phase-space functions. 
This map is also a Lie algebra homomorphism and this is a key step in the following discussion, 
which explains how the operation of taking the moments is a momentum map arising from a 
Lie algebra action on the Vlasov distribution. Indeed we shall see that the momentum map is 
exactly the dual of this Lie algebra homomorphism. 

Analogous considerations [28] on the structure of the moment Lie algebra also hold for the 
Lie algebra of statistical moments, which are more simply defined as symmetric tensors, rather 
than tensor fields. 

3 A momentum map for kinetic moments 

As a starting point, one recalls [33] the definition of a momentum map: Given a Poisson manifold 
(a manifold P with a Poisson bracket {•,•} defined on the smooth functions J-{P)) and a Lie 
group G acting on it by Poisson maps (i.e. the Poisson bracket is preserved), a momentum map 
is defined as a map J : P — > £* such that 



where £ is the Lie algebra of G, the notation (•, •) indicates the pairing between £* and £, and 
£p is the vector field given by the infinitesimal generator 



Thus, whenever the infinitesimal generator is known, a corresponding momentum map can be 
defined. Indeed, the Lie group action is not strictly necessary to this definition. Rather what is 
necessary is the Lie algebra action [44] . i.e. the map associating the infinitesimal generator £p 
to the generic Lie algebra element £. 

In the case of moments, the Lie-Poisson structure (|3|) involves a Lie algebra which is defined 
by the term in square brackets, while the moments A n belong to the dual Lie algebra. Moreover, 
moments are linear functionals of the Vlasov distribution, /, as in ([3]). The aim of this section is 
the show that the operation of taking the moments is not only a Poisson map (Gibbons' theorem 
[20] ). but is also a momentum map, arising from a Lie algebra action on the phase space density 
/. Therefore, the Poisson manifold P in the definition of J is the space of densities P = 
T*(T*Q), while the moment Lie algebra is defined through the expression in (|3|). Consequently, 
a single element £ S £ is identified with a whole sequence (5 n and the corresponding infinitesimal 
generator £p is denoted by (3p. In the remainder of this section the bracket {•, •} always denotes 
the canonical Poisson bracket, while {•, -}y stands for the Vlasov bracket ([2]). 



{F(p),(J(p),0} = Zp[F(p)} 
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3.1 The Lie algebra action and its momentum map 



In order to show that moments constitute a momentum map, one starts with the Vlasov equation 
written in the advective form 

% + £x»/ = 0, 

where the phase-space function H is the single-particle Hamiltonian, which is written from the 
Vlasov Lie-Poisson structure as H = SH/Sf. Substitution of expression ([7]) then yields the 
Vlasov equation in the form 

The implied summation on repeated indices will be omitted from now on. We will also use the 
notation for the one-dimensional case, since the extension to higher dimensions may be madel 
by following the arguments in the previous section. The bracket term is evidently determined 
by the Lie derivative «£x / along the Hamiltonian vector field X defined through the fiberwise 
polynomial p n 57i/5A n . That is, the bracket term is determined by the symplectic Lie algebra 
action on its dual. However, by the usual isomorphism between polynomials and symmetric 
tensors, this term may also be considered as arising from the moment Lie algebra action on the 
phase space densities, so that 

P-f = - {f,P n Pn} , (20) 

where (5 denotes here a whole sequence of symmetric contravariant tensor fields (3 := {f3 n } ne ^- 
At this point, one defines the map J : !F*(T*Q) — > g* (g is the Lie algebra of contravariant 
tensor fields, as above) from Vlasov distributions to moments as in (J3]) (or as in (jlip . when 
including the basis). A direct verification shows that this is a momentum map. Namely, 



P n f(q,p)dp Pn(q)dq } dqdp 



jjf(q,p)^,p n Pn(q)^dqdp 



f(q,p),p n Pn(q)}^ydqdp = (3 P [F(f)], 



Sf 

where the last line is justified by integration by parts and by the definition of the moment Lie 
algebra action in fl20J) (recall that P = F*(T*Q) as above). 



3.2 Infinitesimal equivariance 

This section shows that the momentum map J(/) defining the moments is (infinitesimally) 
equivariant |44j . i.e., 

/3 P [(j(/),7>] = (ad£ J(/), 7 > V/?, 7 eg V feP = F*(T*Q), (21) 

where the coadjoint operator ad* is defined as in (I16p . Although one can easily verify this 
property by direct substitution upon following the calculations in [201 [231 [23] , we prefer to show 
equivariance through the identification of J with the dual a* of a Lie algebra homomorphism a, 
which is known [H] to be an equivariant momentum map in the general case. This is preferable, 
since it enriches the geometric character of J. 

A Lie algebra homomorphism is a map a : TO — ► f) between two Lie algebras tv and f) satisfying 
the following property 

Vlto) = H6, «(??)](, • 
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In the present case, one has to = g (the moment algebra) and f) = J-(T*Q). The homomorphism 
a is actually provided by the usual isomorphism between symmetric tensors and polynomials, 
so that a : j3 \— > p n f3 n . The Lie brackets on g and T are defined by the Schouten concomitant 
and the canonical Poisson bracket respectively. A trivial calculation shows that 

«([/?, 7]) = P m+n ~ 1 [A»,7m] = {P n pn, V m lm} = {«(/?), «( 7 )} ■ 

Consequently, 

<J(/),/3) fl = jA n (f)P n dq = jjfp n (3 n dqdp=(f, a (f3)) :F , 

so that J = a* . The (infinitesimal) equivariance (|2ip of J arises naturally, since it is a general 
fact [44] that any dual Lie algebra homomorphism is an equivariant momentum map. Notice 
that from (|2ip one can characterize equivariance also by 

«J(/),/3), <J(/), 7 >} V = <J(/), [P, 7]), 

which is verified by direct substitution. At this point, Gibbons' theorem on Vlasov kinetic 
moments becomes justified by the equivariance of the momentum map J (/) = {^4n} n 6 n ■ hideed, 
any equivariant momentum map is Poisson [44] . 

The term "infinitesimal" used above refers to the fact that the momentum map is not known 
to arise from a group action. Indeed, the group action associated to moment dynamics is still 
unknown and we leave this issue for future work. In particular, one would like to identify 
the moment algebra q as the tangent space at the identity T e G of some Lie group G, whose 
infinitesimal action is given by (|20p , 

4 A bosonic Fock space for the moments 
4.1 Bosonic structure of the moment Fock space 

It is important to notice that if one fixes a point q on the base manifold Q, then the /3's can 
be considered as simple tensors belonging to V T q Q, the symmetric algebra on the tangent fiber 
T q Q. Here the notation \f W stands for (§) n (\/ n W), where V denotes the symmetric tensor 
product |51| . so that \J n W := S ((^) n W) and S is the symmetrizing operator. Thus, the space 
of moments can be identified with the Fock space on the tangent fiber T q Q corresponding to the 
Bose-Einstein statistics |17j . 

Remark 2 (The bundle of Fock spaces) In more generality, the moment algebra is com- 
posed of tensor fields, mapping the configuration manifold Q to the bundle of Fock spaces 

STQ := (J (V^Q) • 

The geometric properties of such bundle are not completely understood. An important question 
would concern its dual ST*Q := U g (\J T*Q) . Indeed, it is well known that the cotangent bundle 
T*Q = [J q T*Q is characterized by the symplectic form lo = d8 = d(pdq) and it is not clear what 
role (if any) this symplectic form may play on ST*Q. Such questions will be addressed in future 
work. 
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4.2 Creation and annihilation operators 

In second quantisation one also defines creation and annihilation operators. Within the moment 
Lie algebra, one may choose a vector field a\ and use it to construct the creation operator acting 
on a symmetric n-tensor field (3 n as 

a(ai)(3 n :=aiV [3 n := S( ai ®f3 n ) . (22) 

The annihilation operator can be defined as the dual 

(A n+1 , a(ai)(3 n ) = ((ati-lSA n+1 ), (3 n ) =: (a*(ai)A n+l , (3 n ) , (23) 

so that 

a*(a 1 )A n = (a 1 -lA n ) , (24) 

where the contraction is consistently symmetrized. Consequently these operators are such that 
a:T q Qx V™ T q Q \/ n+1 T q Q and a* : T q Q x V" T*Q -» V™" 1 T*Q. It is obvious that these 
operators commute 

[a( ai ), 0(71)] = = [a*(ai), a*( 7l )] . (25) 

The fact that a acts on symmetric powers of the tangent space is not essential: one only needs to 
think of a and a* as operations a : W x \J n W -»■ Wand a*:lf x V™ ^* -» V"" 1 ^ EH 

whose properties are unchanged when W can be identified W*. In particular, if one considers 
the moment tensor-densities A n of the form A n = Tl n <8> Aq (with LTo := 1), one can construct the 
operation 

a(0i)A. ==©l V^n := (6i vn„)®i (26) 
for any one-form Gi € T*Q (one has a(Bi)Ao := 9i <8) Ao). By analogy, one writes the dual 

a*(ei)/3 n = (0i_l/3 n ) . (27) 

By using these definitions one checks that 

[a* (ai), o(0i)] = (ax_iei)I. 

If one takes a\ = dh and Q\ = dq k , then 

[ a hi a k] = $hk ! 

where one defines a& := a (d<7 fc ) and aj^ := a*((?/ l ). These commutation relations are identical to 
those used in second quantization. 

It is important to notice that any one form Gi can be identified with a first-order kinetic 
moment B\ = @± (8) -Bo arising from a Vlasov distribution. In more generality, any element 
G n G V™ TqQ can be identified with a n-th order moment B n = Q n <g> Bo, so that the operators 
a(@i)A n and a*(@i)Pi can be considered as operators such that 

aiB^An := ( |i V A n ) ® ^ o*^)/^ = ^ _l A, 

\Bq J iJo 

At this point, creation and annihilation operators involve only moments A n ,B n and their dual 
variables a n , /? n . 

It may be useful to point out that the physics literature on quantum mechanics [TTj uses 
the opposite notation for the creation and annihilation operators, so that they become switched 
a <-> a*. In this paper we adopt the more mathematical convention in [51] and no confusion 
should arise from this choice. 
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4.3 Cold plasma solution from creation operators 

A particular case where the creation operator can be applied is the cold plasma solution of the 
Vlasov equation, which is written as f(q,p,t) = Ao(q,t) 5{p — U(q,t)). This yields a moment 
hierarchy such that 

It is now clear that this moment hierarchy may be obtained by iterating the composition of the 
creation operator 

A n = [a(A 1 )] n A . 

Thus, this case provides a useful interpretation of creation and annihilation operators in the 
moment Fock space. Any moment is constructed by recursively applying the creation operator 
on the vacuum state Aq, so that A± = a(^4i)^4o = II <g> Aq. 
By analogy, for a vector field ol\ we can see that 

a*{ ai )A n = Oi J IT 1 ) ® A = (ai _l n) A n . x 

and, thus, 

where the sharp Jj indicates the usual operation of raising indexes and |I1| 2 := n_Hl". 

4.4 Occupation numbers 

The analogy with the bosonic Fock space in second quantisation proceeds further by noticing that 
the multi-index notation introduced in equation ([5]) corresponds to the treatment of occupation 
numbers o"i, . . . , cr/v, f° r the symmetric tensors A n {q) (cf. e.g. [51]). Thus, the interpretation 
in terms of symmetric tensor fields is also equivalent to the geometric characterization of the 
moments A a . The operators a and a* also have an equivalent in terms of occupation numbers. 
For example, take a Lie algebra variable (3 a = (3 ai ,...,a N with ^ cij = n and a vector field ol\,. 
The creation operator becomes 

a (aij Pa = («i V Pn) au „^ k+lr __ )aN , 

so that the £:-th occupation number <7fc is increased by 1. In particular, one may focus on the 
basis and introduce the ket notation 

\ai, . . .,a N ) := {d q i) ai . . . {d q N) aN = \a{) . . . \a N ) , 

so that 

ah |ci, • • • ,o~n) = o(c^) | (Ti, • • • ,crjv) = y/a h + l \a±, . . . ,a h + 1, . . . ,a N ) , (28) 
a% \ax, . . .,a N ) = a*(dq k ) \a 1: . . .,a N ) = y/a^\ai, . . . ,a k - 1, . . .,a N ) . (29) 

Indeed, it is well known in quantum theory [16] that all the quantum- mechanical properties 
of these operators follow uniquely from the commutation relations in the previous section. By 
proceeding analogously, one defines the bra corresponding to the moment variable as 

(<ti,...,<t N \ ■= {dq 1 )* 1 ...(dq N ) aN ®A dVol. 



A n = \U\ 2n A , 
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(The "vacuum density" Aq is not a problem since it yields 1 under the L? pairing.) One then 
writes 

(<7i, ...,a N \ a k = (ai, ...,a N \ a(dq k ) = \Jo k + 1 \a\, ...,a k + l,.. .,a N ) , (30) 

(ai, ... ,<tjv| a* h = (<r\, ... , crjvl a*(d^) = y/ak\ai, . . . ,a h - 1, . . . ,a^) . (31) 

Thus, the moment framework on the bosonic Fock space is completely consistent with the 
standard treatment of second quantization. For example, one can also define the number operator 
M := J2k a k a h so that 

M\a h ) = x /aJ l \a h - 1) = a h \a h ) (32) 

k 

and more particularly, 

N 

N\ai, . . .,<tn) = ^ Q"fe = n (33) 

k=l 

so that the operator N corresponds to the moment order, instead of the number of particles 
occupying a certain quantum system. 

Remark 3 (Statistical moments) Since the treatment of statistical moments is analogous to 
that of kinetic moments and it involves the symmetric tensor algebra in a similar fashion, we 
expect no obstacles in transferring the same Fock space treatment to statistical moments. 



5 Extension to anisotropic interactions 

The Vlasov kinetic equation (pQ) regulates multi-particle dynamics in phase space of pointdike 
particles. What happens if the particles have a shape in physical space (e.g., they are rigid bod- 
ies)? This section addresses the question of how to generalize the geometric moment framework 
to particles undergoing anisotropic interactions. We base our analysis on a particular kind of 
Vlasov equation, which was proposed by Gibbons, Holm and Kupershmidt (GHK) [21|. [22] for 
obtaining the equations of chromohydrodynamics of quark-gluon plasmas. The GHK approach 
defines a density variable which depends not only on phase space coordinates, but also on an 
extra degree of freedom g (e.g., orientation). In order to respect the geometric structure of the 
Vlasov equation, this variable is supposed to belong to a (dual) finite-dimensional Lie algebra 
h, so that the GHK- Vlasov equation is written as 



df dH 
dg' dg 



. (34) 



At this point, one may introduce two possible types of moments: the Smoluchowski moments 
(placed in this framework by Holm, Putkaradze and Tronci [29} [30] ) 

A n (q,g,t) := J p n fd N p 

and the GHK-moments [2H E2] 

AnA^t) ■= I p n 9 u fd N pd K g, (35) 
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where n, v S N and K = dimf). Thus g v is again a tensor z^-th power of g. The approach in 
the previous sections can be generalized to both of these kinds of moments. However, since the 
GHK-moments involve also powers of g they appear to be geometrically more interesting and 
we choose to focus on them for the remainder of this section. 

In order to generalize the Schouten bracket, one starts with the Lie-Poisson structure 



{F, G} [/] = 
The chain rule allows one to write 



^ ^7 '^7 ^ 



d_6F_ d_6G 
dg Sf ' dg 5f 



dq dp d K g 



(36) 



(37) 



and a direct substitution into equation (J36J) yields 



{F,G} 



A 



m.+n— 1, u+fj. i 



imi fin,. 



+ A 



m+n, 



a, 



C 



where [•, •] is the ordinary Schouten bracket, while 



\C 



is the Lie bracket defined by 



a u , I3a] c := u t jL C - 1 ( a » ® Pp) 



(38) 



(39) 



where ® denotes the ordinary tensor product (not symmetrized) and C is the structure tensor 
of f) with components C^ c . The whole bracket tensor [-,-]c 1S then identified with its symmetric 
part by the usual argument. The derivation of (f38|) is presented in appendix 1X1 

Thus, the moment dynamics is extended to anisotropic interactions. It is straightforward 
to verify that truncating the hierarchy to consider only Aq q, Aq \,A\ q yields the equations of 
chromohydrodynamics in [211 122j . Moreover, it is easy to see that the geometric considerations in 
the previous sections can be suitably generalized to the anisotropic case. In particular, moments 
are now defined as tensors 



-4, 



'il,. 



! air 



dq 11 <g> 



dq 1 



3 oi 



d 3 q, 



where {e 1 , . . . ,e A } is a basis of {)* . Thus these moments are defined on the product of Fock 
spaces \JT*Q (g) and all the second quantisation techniques can be as well applied to this 
case. A natural case for future study could be a system of rigid-body particles. In this case 
one identifies a = so (3) ~ K 3 and the structure constants are given by the Levi-Civita symbol 
Eijk corresponding to the vector product in IR 3 . We leave the geometric investigation of this 
structure as a possibility for future work. 



6 Conclusions and open questions 

In this paper, the Schouten symmetric bracket has emerged for the first time as a basic object 
for the applications of Vlasov kinetic theory. In its previous treatments, the Schouten bracket 
has been appreciated as an interesting invariant differential operator, irrespective of its potential 
applications in physical problems. However, many potential applications for it may now arise 
in phase space dynamics, whose primary physical examples are kinetic equations. The key to 
seeing the importance of the Schouten bracket for these applications is the identification of 
kinetic Vlasov moments with tensor densities, rather than simple covariant tensors as in the 
previous physical literature [25] . 
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Directions for further investigation of the Schouten bracket naturally arise from the results 
in this paper. For example, the identification of the moments with the bosonic Fock space shows 
how the cold plasma closure is constructed using standard methods in second quantization, in 
which the vacuum state is given by the particle density Aq. In more generality, a bundle of 
Fock spaces can be defined in order to treat moments as tensor fields, rather than tensors. The 
investigation of the geometry of this bundle is another open question that deserves to further 
investigation. The most remarkable result of the Fock space treatment is the interpretation that 
it reveals of Kupershmidt's multi-indices as occupation numbers. This natural interpretation 
unifies the two approaches and places them into a single geometric framework. 

The present geometric treatment also characterizes the moment hierarchy as a momentum 
map. This property was not entirely unexpected; but it is welcome, because the Lie algebra 
action (|20f) can now be uniquely defined as an action of symmetric contravariant tensor fields, 
so that the momentum map is a dual Lie algebra homomorphism. In this setting, it is clear that 
the Poisson structure for the moments arises from infinitesimal equivariance of the momentum 
map. The outstanding question that remains concerns the characterization of the Lie group 
corresponding to the moment Lie algebra. If this were possible, one could characterize moment 
dynamics as a form of coadjoint motion [H] 



where Ad* indicates the coadjoint group action and A denotes the entire sequence A := {A n } n ^. 
This point is related to the question of global equivariance |H] and would be crucial for the 
identification of families of moment invariants analogous to the Poincare invariants of the particle 
phase space. This identification is known for the statistical moments [28J, but it remains to be 
discovered for the kinetic moments. 

Finally, the geometric framework for the moments has been extended here to apply in the 
case of anisotropic particle interactions. Interest in this extension arose from the previous use 
of kinetic moments in problems of quark-gluon plasma dynamics |21[ [22] . Additional interest 
has also arisen recently in physical problems involving the self-assembly of oriented nano-rods 
[29[ 130] . This extension of moments to the kinetic theory of particles undergoing anisotropic 
interactions may also be useful in problems involving the dynamics of complex fluids possessing 
anisotropic order parameters |26| . 
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A Evaluation of the moment brackets 

This appendix presents the evaluation of the moment Lie- Poisson bracket in the anisotropic case 
(|38p . This proceeds by a systematic full calculation involving tensor indexes [53] • The first term 
in (I38j) is identical to the ordinary moment bracket (113j) . which is thus also evaluated in this 
appendix. One takes the expression in ([36]) and substitutes ([37]) with 5F/5A mjU = a m> „ and 
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5G/5A nifl = P n>fl . At this point one follows the steps below 

{G,H}= Ijjf {<w(q) _| p m ® g", & jM (q) _l p" ^ d^q A d^p A d*g 



d(a m )^-^d Pjl ... Pjn . jlr .. Jn 



dpk 



.4 



dq h dph 
— (a m ,„(q) _l p m ® V ) , (Pn,M) -I P n ® »' 



m+n— 1, n+v 



■l> 



d N q 



+ HI fp ew(q)/?n !M (q) ( 9dC bc - 



dgb 



dg c 



d^q A d^p A d K g 



d^qAd^pAd^ 



.4 



m+n— 1, /i+f 



■I' 



d 3 q 



- ////P m+n (a m ,,) ai --' a "- 1 ' (/3„, M ) S 

(hk g d Ct g ai ... g av _ Y g Sl . . . d^q A d^p A d^g 



.4 



m+n— 1, /i+f 



d 3 q 



.4 



m+n— 1, u+1/ ) 



C^m,;/! An,/ 



+ ! / |j Um+n,^i/-l , C _l a m ,u ® An,/ 



• ^-^-m+n— l,u+i/ ) Qm,ui Pn,fi \ "I - (^^4 



m+n, 



C^m,ui A 



n,u 



which finally coincides with equation 
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